library(pracma)

## 第1题
# A=matrix(c(2,1,2,4),nrow=2,byrow=T) #按行输入矩阵
# print(A)
# E=diag(2) #单位矩阵
# B=A%*%A-2*A+3*E
# print(B)

##第3题
# A=matrix(c(1,2,1,3,1,4,1,5,1),nrow=3,byrow=T)
# B=matrix(0,nrow=3,ncol=3)
# for (i in 1:3){
#   for (j in 1:3){
#   B[i,j]=(-1)^(i+j)*det(A[-i,-j]) #计算代数余子式
#   mystring=sprintf('A%d%d = %3.0f',i,j,B[i,j])
#   print(mystring)
#   }
# }
# round(A%*%t(B),1) # 验证

##第4题
# A=matrix(c(1,2,0,3,0,1,5,11,0),nrow=3,byrow=T)
# E=diag(3) #创建三阶单位矩阵
# AE=cbind(A,E) #左右连接两个矩阵
# print(AE)
# print(rref(AE)) #计算行最简形

##第5题
A <- matrix(c(3, 1, 2, 1), nrow = 2, byrow = TRUE)
B=matrix(c(2,5,5,2),nrow=2,byrow=T)
AB=cbind(A,B) #左右连接两个矩阵
print(AB)
print(rref(AB)) #计算行最简形


